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Abstract 

We investigate tlie phase diagram and, in particular, the nature of the the 
multicritical point in three-dimensional frustrated A^-component spin models 
with noncollinear order in the presence of an external field, for instance easy- 
axis stacked triangular antiferromagnets in the presence of a magnetic field 
along the easy axis. For this purpose we study the renormalization-group flow 
in a Landau-Ginzburg- Wilson 0^ theory with symmetry 0(2)®[Z2®0(iV — 1)] 
that is expected to describe the multicritical behavior. We compute its MS 
(3 functions to five loops. For > 4, their analysis does not support the 
hypothesis of an effective enlargement of the symmetry at the multicritical 
point, from 0(2)®[Z2®0(iV - 1)] to 0(2)®0(iV). For the physically interest- 
ing case A'^ = 3, the analysis does not allow us to exclude the corresponding 
symmetry enlargement controlled by the 0(2)00(3) fixed point. Moreover, 
it does not provide evidence for any other stable fixed point. Thus, on the 
basis of our field-theoretical results, the transition at the multicritical point 
is expected to be either continuous and controlled by the 0(2) 00 (3) fixed 
point or to be of first order. 

PACS Numbers: 64.60.Kw, 05.10.Cc, 05.70.Jk 
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I. INTRODUCTION 



Frustrated spin models with noncoUincar order, such as easy-axis stacked triangular 
antiferromagnets (STAs), reveal a quite complex phase diagram in the presence of an external 
magnetic field H along the easy axis, see, e.g., Refs. [1-9]. A model of these systems is 
obtained by considering a stacked triangular lattice, three-component spins Sj defined at the 
sites of the lattice satisfying Sj ■ Sj = 1, and the Hamiltonian 
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JijSi ■ Sj + Y^[Dsl^ + Hsi 



[1.1) 



with an antiferromagnetic hopping term Jij. For small magnetic fields, one observes two 
critical lines that are expected to belong to the XY universality class according to the 
theoretical analysis [10-13]. For large magnetic fields, there is instead a single critical line 
that is expected to belong to the 0(2)(S)0(2) universality class, which is characterized by 
the symmetry-breaking pattern 0(2)®0(2) — > 0(2)diag- Finally, the large-if and small-if 
domain are separated by a first-order spin-fiop fine. These four critical fines meet at a 
tetracritical point, see Fig. 1. The authors of Ref. [12] argued that the critical behavior at 
the multicritical point belongs to the 0(2)(S)0(3) universality class, with symmetry- breaking 
pattern 0(2) 00 (3) — >• Z,2®0(2)diag, described by the Hamiltonian [14] 
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where is a 3 x 2 matrix, i.e., a = 1, 2, 3 and i = 1,2. The analysis of Ref. [12] is however 
not complete, since only the quadratic perturbations of the 0(2) (8)0 (3) fixed point (FP) 
were considered. In this paper we reconsider the issue, by performing a complete analysis of 
all quadratic and quartic perturbations induced by the easy-axis anisotropy, including also 
those terms that arc absent in the theoretical analysis of Ref. [12]. 

The Landau-Ginzburg- Wilson (LGW) Hamiltonian that describes the multicritical be- 
havior is the most general Hamiltonian with symmetry 0(2)®[Z2©0(2)]. It is given by 
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where (pai and ipi are real fields with a = 1, 2 and i = 1,2. The LGW Hamiltonian (1.3) can 
also be obtained from model (1.1) by performing a Hubbard-Stratonovitch transformation 
[15,16] and an expansion in terms of the critical modes. 
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FIG. 1. The experimentally observed phase diagram in easy-axis STAs as a function of tem- 
perature T and magnetic field H. 

The critical behavior at the multicritical point is determined by the stable FP of the 
renormalizat ion-group (RG) flow when both and are tuned to their critical value. 
See, e.g., Refs. [17,18] and references therein for a discussion of the field-theoretical (FT) 
approach to multicritical phenomena. If no stable FP exists or if the system is not in the 
attraction domain of the stable FP, the transition at the multicritical point is expected to 
be of first order. The hypothesis of the effective enlargement of the symmetry 

0(2) ® [Z2 © 0(2)] ^ 0(2)® 0(3) (1.4) 

at the multicritical point requires that the 0(2)®0(3) chiral FP is stable with respect to 
the quartic terms that break 0(2)®0(3) to 0(2)(8>[Z2©0(2)] — those proportional to Wq, yo, 
and Zq in Hamiltonian (1.3). If this does not occur, the 0(2)®0(3) FP does not control the 
multicritical behavior for generic values of the Hamiltonian parameters. As a consequence, 
the effective enlargement of the symmetry to 0(2)(8)0(3) at the multicritical point requires an 
additional tuning of the parameters: beside tuning and r^, at least one more Hamiltonian 
parameter must be properly fixed to decouple the additional relevant interaction. 

In this paper we investigate this issue by FT methods. We consider the more general 
theory in which the order parameter (pai is an {N — 1) x 2 matrix, i.e., a — 1, ...N — 1 
and i — 1,2, for N > 3. In this case, setting wq = yo = Zq = and = r^, one 
recovers the 0(2)©0(iV)-symmetric LGW Hamiltonian (1.2). This theory has a stable FP 
with attraction domain in the region (72,0 > describing a critical behavior with symmetry- 
breaking pattern 0(2)(g)0(iV) 0(2)(g)0(iV - 2) [14,19,20]. We should mention that the 
existence of this FP has been a controversial issue for quite a long time; the different scenarios 
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are reviewed in Rcfs. [2,20-23]. 

In order to determine the RG flow of the theory (1.3) in three dimensions and deter- 
mine its multicritical behavior, we consider the minimal-subtraction (MS) scheme without 
e expansion (henceforth indicated as 3d-MS scheme) in which no e expansion is performed 
and e is set to the physical value e — 1 [24] . We use a symbolic manipulation program that 
generates the diagrams and computes symmetry and group factors, and the compilation of 
Feynman integrals of Ref. [25]. This allows us to compute the 3(i-MS /^-functions to flve 
loops in the full model. The perturbative series are used to determine the FP structure 
and, in particular, to investigate the existence of stable FPs that may describe the critical 
behavior at the multicritical point. The 3d-MS scheme is particularly convenient for the 
three-dimensional FT study of the multicritical behavior. Indeed, the multicritical theory is 
simply obtained by setting = = 0, i.e. by considering the massless theory. Note that 
this is not correct in the three-dimensional massive zero- momentum (MZM) scheme [26] , in 
which a proper tuning of and is needed. 

We summarize the main results of this paper. In order to check the hypothesis of the 
effective enlargement of the symmetry 

0(2)®[Z2©0(iV-l)] ^ 0(2)®0(iV) (1.5) 

for generic A^-component systems at the multicritical point, we study the stability properties 
of the 0{2)^0{N) chiral FP with respect to all quadratic and quartic perturbations that 
are symmetric under the reduced symmetry 0(2)®[Z2©0(A'^ — 1)]. The analysis of the 
corresponding five-loop series does not support the stability of the 0{2)^0{N) chiral FP 
for any value > 4, with increasing confidence as N increases. For N — 3 the results are 
not conclusive. Our FT results do not allow us to establish the stability properties of the 
0(2) (8)0 (3) FP, which may be either stable or unstable. In the former case, the multicritical 
behavior would be controlled by the 0(2)®0(3) FP if the transition is continuous. In the 
latter case, we note that the crossover exponent should be very small, ^4^4 < 0.1. Therefore, if 
the effective quartic Hamiltonian parameters that break the 0(2)(8)0(3) symmetry are small, 
the crossover from the preasymptotic 0(2)00(3) critical behavior to the eventual asymptotic 
behavior is expected to be very slow, and one may observe an effective 0(2) (8)0 (3) critical 
behavior for a wide range of reduced-temperature values. We also perform a general Pade- 
Borel analysis of the RG flow for = 3, to investigate the existence of stable FPs for generic 
values of the quartic couphngs. No evidence of additional stable FPs is obtained. Therefore, 
according to our FT results, for N = 3 the multicritical transition is either controlled by 
the 0(2)(g)0(3) FP or is of flrst order. 

The paper is organized as follows. In Sec. II we show how the general Hamiltonian 
(1.3), that has been written down on the basis of symmetry considerations, can be recovered 
from Hamiltonian (1.1) of easy- axis STAs in a magnetic held along the easy axis. Then, 
we discuss the mean-fleld phase diagram, showing that, beside the tetracritical behavior 
that was predicted in Ref. [10], the model also admits a bicritical and a pentacritical phase 
diagram. In Sec. Ill we focus on some particular cases. In Sec. Ill B we discuss the stability of 
theO(2)®0(A^) FPs, and, in Sec. IHG, the stability of the decoupled [0(2)(g)0(Ar-l)]©0(2) 
FPs. In Sec. IV we discuss the full model: in Sec. IV A we consider the one-loop e expansion, 
while in Sec. IV B we numerically investigate the RG flow for N — 3. In Sec. V we present 
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our conclusions and critically discuss the experimental results in view of our findings. In 
App. A we report a discussion of the mean-field diagram, in App. B we classify all quadratic 
and quartic perturbations of the 0(2)(8)0(iV) FPs that are 0(2)-invariant. Finally, in App. C 
we compute the RG dimensions of all quadratic perturbations of the 0(2)®0(A^) symmetric 
theory. 



II. DERIVATION OF THE GENERAL HAMILTONIAN AND MEAN-FIELD 

ANALYSIS 



In this Section we derive the effective LGW Hamiltonian (1.3) for easy- axis STAs in a 
magnetic field along the easy axis described by Hamiltonian (1.1). It can be obtained as 
usual by first performing a Hubbard-Stratonovitch transformation [15,16]. If a = 1, 2, 3, 
is an unconstrained three-component field, the partition function can be rewritten as 
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Jd^s5{s'^ - 1) exp{-Dsl - Hs,) 
As usual, we now expand iy($) in powers of At order we can write the expansion as 



W{^) = an$^ + a2i{^l + $J) + 022$' + 031$^ + a^2^,{^l + 



(2.3) 



where be easily computed in terms of H and D. Finally, we identify the critical 

modes that are associated with the wavevectors Q that maximize the kinetic term. In our 
case the relevant modes are associated with the wavevector Q — (47r/3, 0, tt), so that we can 
write 



$x(r) = c[0ii(r) cos(Q ■ r) + 0i2(r) sin(Q • r)], 
^y{r) = c[(j)2i{r) cos{Q ■ r) + (j)22{r) sm{Q ■ r)], 
^zir) = c['ipi{r) cos{Q ■ r) + il)2{r) sin((5 ■ r)]. 



(2.4) 



where we have introduced new fields (pai and ipi^ with a = 1,2 and i ~ 1,2. The constant c 
is fixed so that, keeping only slowly varying contributions, we have 



ij \_ ai i 



■ (2.5) 



Substituting expressions (2.4) in Eq. (2.3) and keeping only slowly varying contributions, 
we obtain the Hamiltonian (1.3) with = (c^a2i + ao), = (0^022 + ^o); = 9a42c''/2, 
vq — 3a42C^, wq — 9(2a4i — 042)0^/2, yo = 9(2043 — a42)c''/2, and zq — 3(2a4i — 042)0^. This 
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derivation gives relations among the different couplings appearing in Eq. (1.3). However, 
they should not be taken seriously, since the effective Hamiltonian (1.3) is only an approx- 
imation of the original one. Note that the odd powers of $ present in Eq. (2.3) do not 
contribute to the effective Hamiltonian. Indeed, here the basic ingredient is the anisotropy, 
i.e. the breaking of the 0(3) spin symmetry to Z2©0(2). The additional breaking of the 
Z2 symmetry caused by the magnetic field does not play any role, apart from modifying the 
explicit expressions of the parameters of the LGW Hamiltonian. 

Several comments should be made on the derivation of the effective Hamiltonian (1.3). 
First, we have assumed here that the relevant modes can be inferred from the analysis of the 
hopping term. While this is correct for unfrustrated systems, for frustrated ones as is the 
case here the method is questionable and it is indeed possible that some low-temperature 
properties are not correctly described by this approach [27]. Second, note that the explicit 
value of Q does not play any role. This means that the effective Hamiltonian (1.3) can be 
used to describe any critical behavior with a three-component order parameter S'(g), with 
q 7^ K/2, where X is a reciprocal-lattice vector. The effective Hamiltonian (1.3) differs from 
that presented in Refs. [10,12], in which the quartic breaking terms proportional to wq, yo, 
and Zo are absent.^ In any case, even if absent in the microscopic model, these additional 
terms would be generated by RG transformations. 

The phase diagram of the theory with Hamiltonian (1.3) can be studied within the mean- 
field approximation. The general discussion is presented in App. A. Here we only report the 
final results for the specific case in which one of the transitions is a chiral transition, i.e. it 
is associated with the symmetry-breaking pattern 0(2)(S)0(2)-^0(2)diag, as it is of interest 
for easy-axis materials. This occurs when vq + zq > 0. The reader interested in systems with 
a coUinear/paramagnetic transition is referred to App. A. The possible phase diagrams are 
reported in Fig. 2. There are three possibihties: (a) a pentacritical point, (b) a tetracritical 
point, (c) a bicritical point. In all cases but one the transitions are second-order ones; one 
transition line is of first order. The known easy- axis materials, like ANiXs, with A = Cs, 
Rb, and X = CI, Br, or CsMnIs, all show a tetracritical point, i.e. a phase diagram of type 
(b). This fact should be related to the smallness of the easy-axis anisotropy. Indeed, for 
D ^ we also have if — > at the multicritical point. Thus, the breaking of the 0(3) 
invariance is expected to be small at the multicritical point, and we can generically assume 
that Wo ~ yo ~ -2o ~ 0. In this specific case, the mean-field analysis predicts a tetracritical 
phase diagram irrespective of uq and vq, in agreement with experiments. Bicritical behavior 
is expected for \D\ large enough: for \D\ > 3J', where J' is the intraplane coupling, the 
basal spin components and Sy should not magnetize at any temperature [1] and thus 
phase 4 should not occur, forbidding a tetracritical behavior. For H = 0, these systems 
should behave as Ising antiferromagnets. Note that, in this case, beside the XY transition 
predicted by the phase diagram (c), other transitions (probably first-order ones) may occur 
as T is lowered [28-30]. It is interesting to note that the mean- field analysis also predicts 



^Hamiltonian (1.3) is recovered by adding in the Ansatz of Ref. [10] for the free energy, see their 
Eq. (3), terms proportional to \S^\^, S'^\S^\'^, and {S^S* ■ S* + {S*)'^S-S). These terms are a priori 
expected in the presence of easy-axis anisotropy on the basis of symmetry considerations. 
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FIG. 2. The possible phase diagrams in the {r^,r^) plane predicted by the mean-field approx- 
imation with a chiral ("ch") transition. Thin lines represent second-order transitions, while thick 
lines are first-order transitions. Phase 1 is paramagnetic, in phase 2 ijj ^ {) and = 0, in phase 3 
V' = 0, </>i / 0, 02 / with 01 • (^2 = 0, in phase 4 0i 7^ 0, 02 = 0, 7^ 0, in phase 5 all vectors are 
nonvanishing. 

a pentacritical point. In this case there is a new phase (phase 5 in Fig. 2) in which the 
basal spin components and Sy show a distorted 120° structure, while the component is 
modulated as in phase 2. The possibility of a tetracritical phase diagram was already noted 
in Ref. [10]. However, in their work, since they only considered the case Wq = = Zq = 0, 
this was the only possible phase diagram. Indeed, the pentacritical and the bicritical points 
are only obtained if Wq, |/o, and Zq are not all vanishing. 

III. ANALYSIS OF SOME PARTICULAR CASES 

A. Particular models and fixed points 

The three-dimensional properties of the RG flow are determined by its FPs. Some of them 
can be identified by considering particular cases in which some of the quartic parameters 
vanish. For example, we can easily recognize: 

(a) the O(i^) -symmetric model is recovered by setting = Uq = vq = wq = 2:0 = {K ~ 
2), r^p = uq = vo = yo = zo = {K ^ 2N - 2), = and vq = wo = yo = zo ^ 
{K — 2N). Results for these theories are reviewed, e.g., in Ref. [22]. 
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(b) the 0(2)(8)0(ir) model, cf. Eq. (1.2) with ^ai being a K x 2 matrix, for = 
and Wo = Ho = Zq = {K = N), and for = and Uq = Vq = Uq = {K = 
N — 1). The properties of these models are reviewed in Refs. [2,20,22,23,31]. In three 
dimensions perturbative calculations within the MZM scheme [19,32] and within the 
3d-MS scheme [20] indicate the presence of a stable FP with attraction domain in 
the region g2fl > for all values of K, except possibly K = Q. For K = 2, these 
conclusions have been recently confirmed by a Monte Carlo calculation [20]. On the 
other hand, near four dimensions, a stable FP is found only for large values of K, i.e., 
K > K^^ 21.80 - 23.43e + T.OOe^ + 0{€^) [14,33-35]. A stable FP with attraction 
domain in the region g2fl < exists ior K — 2 (it belongs to the XY universality 
class) [14], for K — 3 (Ref. [31]), and, as we shall discuss below, for K = 4. Note that 
nonperturbative approximate RG calculations have so far found no evidence of stable 
FPs for i^T = 2 and 3 [23,36]. In the following we will call the FP with g2 > chiral 
FP, while the FP with (72 < will be named collinear FP. 

(c) the 0{2)®0{K) model with K = 2N — 2 for vq = zq = 0. This theory describes the 
multicritical behavior of a model with two order parameters that is symmetric under 
the group 0(2)©0(i^r) [17]. In the case we are interested in, i.e. for A'" > 3 and 
therefore K > A, the stable FP is the decoupled FP, corresponding to = 0, which 
describes a critical behavior in which the two order parameters and ^jJ are effectively 
uncoupled at the multicritical point [18,37]. 

(d) Decoupled 0(2)(g)0(A^ — 1) and 0(2) models for uq — Vq — 0. The corresponding 
stable FP describes two effectively decoupled critical behaviors in the 0{2)<SiO{N — 1) 
and 0(2) universahty classes. 

(e) For Wo = we obtain a multicritical theory with the larger symmetry group 
[0(2)«)0(Ar- 1)]©0(2). 

The FPs of the above-mentioned particular models are also FPs of the enlarged model (1.3). 
Their stability in the full theory can be checked by computing the RG dimensions of the 
additional terms present in the complete Hamiltonian (1.3). 

The analysis of model (1.3) is simplified by the following symmetry transformation. If 
we transform the fields as 

tpi^^^ijtpj, (t>ai^(t>ai, (3.1) 

3 

and the couplings according to 

Uq^Uq- Vq, Vq -Vo, Wo^Vo + Wq, Uo ^ + Ho, Zq 2vo + Zo, 

(3.2) 

we reobtain the same Hamiltonian. This transformation leaves invariant the model with 
f = and maps any model with < to a model with Vq > 0. The transformation 
(3.2) implies several symmetry properties for the RG functions. In particular, under the 
transformation of the renormalized quartic couplings 
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u^u — v, V ^ —V, w^v + w, y^v + y, z^2v + z, (3-3) 

the RG functions associated witii tiie exponents are unchanged, while the /5-functions trans- 
form covariantly. As a consequence, each FP with nonvanishing quartic couphng v < is 
mapped into another FP with f > and the same stabihty properties. The symmetry (3.2) 
imphes that there exists another 0(2)(g)0(A^) model beside that reported at point (b) above: 

(b') The 0(2)(8)0(A'") model is also obtained for = r^, wq — —vq, yo — —vq, zq — —2vq; 
in this case 511,0 ^ uq - vq, 512,0 = -^'o- 

Because of symmetry (3.3) it is enough to study the RG flow for v > 0. 

The results for models (a), (b), (b'), (c), and (d) allow us to identify four possible FPs that 
are candidates for being stable FPs of the full theory. In the 0(2)®0(iV) model (1.2) there 
may be, depending on the value of A^, two FPs. The chiral FP is located at gi = glchi-^) 
and g2 = 5*2 ch(-^) with 5*2 ch(-^) > exists for any value of N except possibly = 6 
[20]. The co'Uinear FP is located at gi = g^^iiN) and g2 = g^^iiN) with 5';^(^) < 0- Such 
a FP exists for A" = 2 [14] — in this case it is equivalent to the standard XY FP — and for 
A^ = 3 [31]. We have investigated if the collinear FP exists also for larger values of A^, by 
extending the analysis of Ref. [31]. In the 3(i-MS scheme at five loops, we find a stable FP 
only for A^ = 4: 5*^^1(4) = 0.10(7) and 5*^^1(4) = -1.83(10).^ For AT = 5, we find a FP only 
in 1/3 of the cases that are considered and such a percentage decreases as A^ increases. In 
the MZM scheme at six loops a FP is found for A^ = 4, 5, and 6, and disappears for N > 8. 
For A^ = 4 it is located at ^jf,ci(4) = 6.1(4) and ^2,ci(4) = -50(2); for A^ = 5 and A^ = 6 
at = -7.2(5) and gl^^ = -53(3), gl^^ = -8.5(7) and = -56(4), respectively The 
perturbative analysis provides therefore strong evidence for the existence of a collinear FP 
for 2 < A^ < 4. For N > 8 this FP is absent (in agreement with the large- A^ analysis), 
while in the intermediate cases 5 < A" < 7 it is not clear whether the collinear FP really 
exists since the two perturbative schemes give opposite results. We have also verified the 
stability of the collinear FPs within the 0(2)®0(A^) theory (1.2): whenever they exist, they 
are stable. However, for A^ = 4 we have been unable to estimate the stability eigenvalues. 
Indeed, the 3o?-MS scheme predicts complex stability eigenvalues, while in the MZM scheme 
we find complex eigenvalues only in 50% of the cases. 

The collinear FP has g2 < 0. Therefore, if we are only interested in the full model for 
V > 0, we must consider the FP appearing in model (b'). Thus, the above-reported results 
for the 0(2) (8)0 (AT) predict two possible FPs: 

(a) u = ^i%h(^), V = ^2,ch(^), w^y^z^O; 



^In the Sd-MS scheme the couplings are normalized so that = gifiH~'' / with = 
2'^~^7r''/^r(d/2). In the MZM scheme the renormalized couplings gi and g2 arc normalized so that 
9i = gifi/m at tree level, where m is the renormalized mass. For reference, we report the collinear 
FP for A = 3 in these normalizations (Ref. [31]): g^^^^ = -7.0(5) and g^^^^ = -50(2) in the MZM 
scheme, ^^^.j = 0.04(8) and g^^^ = —1.71(9) in the 3(i-MS scheme. The same normalization is used 
for the MS couplings in Sec. IV A. 
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(c) u = v = 0,w = gl^^{N -l),z = ^2*ch(^ -!),!/ = ^xvi 

(d) u^v^O,w^ gl^,{N -l),z^ ^2%,(7V - 1), y = ^^y- 



Here (yf^y ^^e four-point rcnormalizcd coupling of the XY theory (numerical estimates can 
be found in Ref. [38]). Note that FPs (c) and (d) are also FPs of the theory with v — 0. 



In this section we study the stability properties of the two FPs that appear in the 
0(2) (8)0 (A/") model and, therefore, we check the possibility of an enlargement of the symme- 
try at the multicritical point from 0(2)(8)[Z2©0(A^- 1)] to 0{2)(^0{N). For this purpose we 
need to classify the perturbations that break the O(A^) symmetry to 0(A^ — 1) and preserve 
the 0(2) symmetry, according to their transformation properties under the O(A^) group and 
the number of fields. This is done in detail in App. B. The multicritical Hamiltonian (1.3) 
can be rewritten as 



where Tigym is the 0(2)®0(iV)-symmetric Hamiltonian (1.2) obtained by setting wq = yo = 
zq = and = r^, and V^'^''^\ y (4,2,1)^ ^^^^ ^(4,2,2) respectively a spin-2 quadratic 

term, a spin-4 quartic term, and two spin-2 quartic terms. Their explicit expressions can 
be found in App. B. The spin-2 quadratic pertubation F^^-^) = (p'^/N - {N - l)ip^/N is 
always relevant. Its RG dimension at the 0(2)(8)0(A^) FP gives the crossover exponent, i.e. 
4> = ^y2,2, where u is the correlation-length exponent. 

Let us first discuss the chiral FP (a) that has 5'2ch(-^) > 0- The exponent 1/2,2 coincides 
with the exponent defined in App. C, since 

^(2,2) _ ^fj-/0f2. For iV = 3, we obtain 
^2,2 = 1.49(3) in the MZM scheme and 1/2,2 = 1.54(8) in the 3d-MS scheme. 

In order to estimate the RG dimensions 1/4,2,1, ?/4,2,2, and t/4,4 of the quartic perturbations, 
we computed the corresponding five-loop series in the MS scheme and we analyzed them 
within the 3(i-MS scheme. Note that, since the spin-2 quartic operators mix, 1/4,2.1 and 1/4,2,2 
are the eigenvalues of the corresponding RG-dimension matrix. We also estimated the spin- 
4 RG dimension 1/4,4 by computing the corresponding five-loop series in the MZM scheme. 
In this scheme we cannot estimate the RG dimensions of the spin-2 quartic perturbations 
because they also mix with the lower-dimension spin-2 operator V2,2- Such a mixing does 
not occur in the 3(i-MS scheme, since in this case the theory is massless and, therefore, 
operators of different naive dimensions do not mix under renormalization. 

Here we only report the series for 1/4,4, which will be the most relevant for the analysis 
of the stabihty of the 0(2)(g)0(iV) FP, for iV = 3 and 4, and in the MS and MZM schemes. 
In the MS scheme we have 



B. Stability of the 0(2)(g)0(A^) fixed points 



(3.4) 




(3.5) 
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TABLE I. Coefficients Uij of the five-loop MS and MZM expansions of 2/4^4, of. Eqs. (3.5) and 
(3.6), for AT = 3 and 4. 



i,j 


N = 3 


TV 

Mo 


N = 4 


AT = 3 


MZM 


N = 4: 


1,0 


—2 




—2 


-6/7 




-3/4 


0,1 


-1/3 




-1/3 


—2/9 




—2/9 


2,0 


28/9 




61/18 


0.383976 




O.32O0O2 


1,1 


1/9 




-1/6 


0.018812 




—0.030864 


0,2 


—1/36 




1/8 


—0.005487 




0.041152 


3,0 


—11.157290 




—12.122093 


—0.246143 




—0.176781 


2,1 


-0.245370 




0.984242 


-0.016559 




0.015113 


1,2 


— 1.174531 




—2.872196 


—0.072221 




—0.142869 


0,3 


0.824994 




1.509014 


0.080442 




0.140829 


4,0 


62.535697 




72.040868 


0.225243 




0.152346 


3,1 


-6.135544 




-19.807535 


-0.038487 




-0.083599 


2,2 


14.290714 




31.675613 


0.157437 




0.252196 


1,3 


-8.768847 




-16.889435 


-0.138290 




-0.224249 


0,4 


0.999971 




1.877179 


0.008548 




0.019106 


5,0 


-422.2.34940 




-508.920947 


-0.231515 




-0.138376 


4,1 


90.411170 




230.058460 


0.069042 




0.102807 


3,2 


-166.086612 




-359.525801 


-0.268436 




-0.357188 


2,3 


108.042417 




220.680543 


0.271663 




0.388611 


1,4 


-23.261600 




-49.094226 


-0.061312 




-0.100610 


0,5 


2.329711 




4.789714 


0.015950 




0.019262 



with e = 1 in three dimensions. The coefficients aij for = 3, 4 are reported in Table I to 
five loops, i.e. for i + j < 5. The renormalizcd MS couplings gi and (72, corresponding to 
the quartic parameters gi^ and g2,o of the Hamiltonian (1.2), arc normalized as in Ref. [20], 
see footnote 2. In the MZM scheme we have the analogous expansion^ 

y4,4 = 1 + XI ^ijdlgi (3.6) 

ij 

The coefficients aij to five loops are reported in Table I for A" = 3,4. The renormahzed 
MZM couplings gi^2 are normalized as in Ref. [19], i.e. g^ a; Cigifl/m at tree order, where 
Ci = (8 + A^)/(487r) and C2 = 3/(167r). The perturbative series that are not reported here 
are available on request. 

We analyzed the series using the conformal-mapping method and the Pade-Borel method, 
following closely Refs. [40,20], to which we refer for details. The error on the conformal- 
method results takes into account the spread of the results as the parameters a and h are 
varied (cf. Ref. [40] for definitions) and the error due to the uncertainty of the FP location 
(we use the estimates reported in Refs. [19,32,20,31]). In the Pade-Borel analysis we used 
the Pade [4/1] and several values of the parameter h. Again, the error takes into account the 
dependence on b and the uncertainty of the FP location. In the Pade-Borel analyses of ^4^4 
the error also takes into account the difference between the [4/1] and the [3/1] estimates. 
This was not done for the spin-2 RG dimensions since in that case the results obtained by 



^The perturbative MZM series were obtained by using a symbolic manipulation program that 
generated the diagrams, the symmetry, and group factors. Numerical estimates of the Feynman 
integrals were taken from Ref. [39]. 
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TABLE II. Estimates of the RG dimensions 2/4,2,!) 2/4,2,2) and 2/4,4 of the operators V^^''^'^\ 
and V^'^'^^ at the chiral (ch) FP and at the colhnear (cl) FP. Results have been obtained 
in the 3d-MS scheme (MS) and in the three-dimensional massive zero-momentum scheme (MZM). 
For the resummation of the five-loop perturbative series the conformal-mapping (CM) and the 



Pade-Borel (PB) methods have been used. For 


2/4,4 two estimates are reported 


in each case, 


rcss]:)c-cl 


, ivoly from 


the ana 


lysis of (lie 


sorios of //I I 


and f/yi.i. 






FP.AT 


i/4,2,1 


2/4,2,2 




2/4,4 






MS, CM 


MS, PB 


Ms, CM 


MS, PB 


MS, CM 


MS, PB 


MZM, CM 


ch,3 


-0.7(9) 


-1.8(7) 


0.0(7) 


-0.9(6) 


-0.4(5), 0.12(5) 


-0.4(4), 0.07(3) 


-0.3(2), 0.05(2) 


ch,4 


-0.7(5) 


-1.0(3) 


0.3(3) 


-0.2(3) 


-0.1(2), 0.18(4) 


-0.2(2), 0.12(4) 


0.15(12), 0.16(4) 


ch,5 


-0.8(7) 


-0.6(3) 


0.2(2) 


0.0(2) 


0.0(2), 0.22(4) 


0.02(11), 0.18(5) 


0.3(2), 0.24(5) 


ch,6 


-0.6(6) 


-0.5(2) 


0.25(13) 


0.11(15) 


0.1(2), 0.26(4) 


0.13(10), 0.23(5) 




ch,8 


-0.5(5) 


-0.4(2) 


0.22(9) 


0.17(6) 


0.2(2), 0.34(5) 


0.29(8), 0.32(4) 


0.26(2), 0.29(2) 


ch,16 


-0.3(2) 


-0.17(6) 


0.09(2) 


0.08(3) 


0.56(4), 0.570(9) 


0.60(2), 0.58(3) 


0.56(2), 0.54(2) 


ch,oo 














1 




cl,3 










0.3(4), 0.7(2) 


-0.7(1.1), 0.5(2) 


0.5(1.0), 1.2(7) 


cl,4 










-0.2(8), 0.5(2) 


-1(1), 0.35(15) 


0.5(1.1), 0.9(3) 



using the Pade [3/1] did not look reliable. They varied significantly with b and, for N < 5, 
favored complex estimates of 2/4,2,1 and 2/4,2,2- 

The results of the analyses are reported in Table II. We first computed the spin-2 RG 
dimensions 2/4,2,1 and 2/4,2,2- If Yij, j = 1,2, is the anomalous-dimension matrix, 2/4,2,1 and 
2/4.2,2 are the eigenvalues of Y. In order to determine them, we resummed the elements 
Yij and computed the two eigenvalues. The results we report are obtained by averaging the 
eigenvalues over many different choices of approximants. The results have a quite large error, 
so that it is impossible to draw definite conclusions on the relevance of these operators. For 
iV > 5 the results favor 2/4,2,2 > 0, so that one operator would be relevant, while for N — 3 
it seems likely that the spin-2 perturbations are irrelevant. For the spin-4 RG dimension we 
report two estimates for each case. They are obtained from the analysis of the series of 2/4,4 
and 1/2/4,4, respectively. Their difference should allow us to estimate systematic errors in 
the series resummations that are not taken into account by the spread of the approximants, 
which we usually take as an indication of the error. We find 2/4,4 > for all > 4, with 
increasing confidence as increases (for = 4 we mainly rely on the MZM analysis that 
predicts 2/4,4 > both for the direct and for the inverse series). For = 3, the results of the 
analysis of 2/44 and I/2/44 differ substantially (the error of the estimate from I/2/44 may be 

underestimated) and even have opposite signs. This does not allow us to establish whether 

y(4,4) 

is a relevant or an irrelevant perturbation of the chiral FP. 
Let us briefiy discuss the physical picture in the two cases. If the 0(2)(S)0(3) FP is stable, 
it controls the critical behavior of statistical systems in its attraction domain. Setting 
t = {T — Tmc)/Tmc, where Tmc is the critical temperature at the multicritical point, the 
singular part of the free energy can be written as 

r =\t\^''f{A\t\-'^), (3.7) 

where A is the scaling field associated with the anisotropy — in STA's it will be a combination 
of T, D, and H — u a; 0.6 is the correlation-length exponent of the 0(2)(8)0(3) theory [19,20], 
= 2/2,2i^ ~ 0.9 is the crossover exponent, and f{x) is a scaling function. 
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If the 0(2)®0(3) is unstable with respect to the spin-4 quartic perturbation, one should 
also consider the crossover exponent </)4 4 associated with the spin-4 quartic instability. This 
is expected to be quite small. Indeed, our results indicate 2/4^4 < 0.2 for = 3, so that 
04,4 ^0.1. Therefore, if the 0(2)®0(3) FP is unstable and the effective quartic Hamiltonian 
parameters breaking the 0(2)(8)0(3) symmetry are small, the crossover from the preasymp- 
totic 0(2)(g)0(3) critical behavior to the eventual asymptotic behavior is expected to be 
very slow and one may observe an effective 0(2)00(3) critical behavior for a wide range of 
reduced-temperature values. 

Now, let us consider the coUinear FP (b) that has 5'2,ci(^) < for A/" = 3 and N ^ A. 
For the spin-2 operators the analysis indicates that they are irrelevant at the FP, but is 
not precise enough to allow a quantitative determination of 1/4,2,1 and 1/4,2,2- For N = 3, hy 
using the conformal-mapping method we find that approximately 50% of the approximants 
give complex spin-2 RG dimensions (the real part is always negative) and 50% provide real 
negative estimates. Pade-Borel approximants always give real negative estimates. The same 
pattern is observed for N — A. The results for the spin-4 operator are reported in Table II. 
They suggest that the collinear FP is unstable for both = 3 and = 4. 

In conclusion, the perturbative analysis shows that the collinear FP is always unstable, 
while the chiral FP is unstable for A^ > 4. For A^ = 3 the results are not conclusive and the 
0(2)® 0(3) chiral FP may be either stable or unstable. 



C. Stability of the decoupled [0(2)(^0(iV - l)]eO(2) fixed points 

Two other interesting FPs can be investigated by essentially nonperturbative arguments, 
FPs (c) and (d) discussed in Sec. Ill A. In order to check the stability of these FPs, we must 
determine the RG dimensions at the decoupled FPs of the perturbations 



J d'xcjj^'', Pt = J d'x Y^ofTij, 

ij 



(3.8) 



where O^^f is defined in App. C, cf. Eq. (C4) and = ipiipj — ^dijip"^ . Simple RG arguments 
show that the RG dimension of Pe, which is an energy-energy term, is given by 

yE^— + 3 = — + — , (3.9) 

where and are the critical exponents of the 3-dimensional XY universality class 

{a^ = —0.0146(8) and z/^ = 0.67155(26), see Ref. [38]), while and z/^ are those of 
the 3-dimensional theory associated with the 0(2) (8)0 (A'" — 1) FP. Analogously, for the RG 
dimension of Pt we obtain the relation 

yT = yr,v + y4,,/> - 3, (3.10) 

where i/t,^ = 1/t,xy, yr,XY is RG dimension of the operator Tij at the XY FP (|/t,xy ~ 1-77 
[41,42]), and 1/4,^ is the RG dimension of the quadratic operator O^f introduced in App. C, 
computed at the 0(2)®0(A^ - 1) FP. 
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Let us consider first the chiral FP (c) that has z > 0. Estimates of the exponents 
and U0 for several values of N are reported in Refs. [19,32,20], while estimates of y4^^ are 
reported in App. C. It is then easy to check that > for iV < 6 and t/e < for iV > 6. 
Therefore, the perturbation Pe is relevant for N < 6. Analogously, we obtain ~ for 
N — 3 and < for A?" > 3. Thus, except possibly for N — 3, Pt is always irrelevant. 
These rcsTilts indicate that the decoupled FP (c) is unstable for A" < 6, and therefore for 
the interesting case A^ = 3. 

Let us now consider the coUinear FP (d). For A^ = 3 the coUinear FP in the 0(2)®0(2) 
theory belongs to the XY universahty class [14]. It is easy to show that = axY, which 
gives ue ~ —0.04. Thus Pe is irrelevant. As for Pt, using the results of App. C, we have 
y4,4> = 2|//i,xY — 3, where yh,XY = (5 — ?7xy)/2 is the RG dimension of the field in the XY 
model. Therefore, in three dimensions yx = |/t,xy — V^^y — 1 ~ 0.73, which means that Pt 
is relevant, and the coUinear decoupled FP (d) is unstable. For A^ = 4 and A" = 5, using 
the results of App. C, we have yr > in both cases. In order to compute yE we determined 

in the 0{2)®0{N - 1) theory at the coUinear FP. In the 0(2)(8)0(3) case we have [31] 

= 0.63(8) (3(i-MS, conformal mapping) and = 0.59(4) (MZM, conformal mapping); in 
the 0(2)®0(4) case we have z/^ = 0.76(9) (3(i-MS, conformal mapping) and z/^ = 0.64(7) 
(MZM, conformal mapping). These estimates imply > for A" = 4 and —0.3 ^yE % 0.2 
for A^ = 5. 

In conclusion, the coUinear decoupled FP is always unstable, while the chiral decoupled 
FP is unstable for N <Q, stable in the opposite case. 

As we have mentioned in Sec. Ill A the decoupled FPs are also FPs for the multicritical 
theory with vq = 0. The analysis of their stability in this particular case follows from the 
results reported above. Indeed, for the theory with vq = it is enough to consider the 
perturbation Pe- Thus, the decoupled chiral FP (c) is unstable for A^ < 6 and stable in the 
opposite case. The decoupled coUinear FP (d) is stable for A" = 3 and possibly for A" = 5. 



IV. THE RENORMALIZATION- GROUP FLOW IN THE FULL THEORY 

A. Renormalization-group flow near four dimensions 

In Sec. Ill we considered some particular FPs, checking whether they were stable in the 
full theory. In order to investigate the presence of other FPs, we now perform a one-loop 
e-expansion analysis. The one-loop /? functions in the MS scheme are: 

N + A 2 N-1 f v\ [N N-l \ 

+?/( 3«- g ) ^(«-^^), (4-1) 

N-6 n ^ wv yv z(N-2) ,^ , 

6 3 3 d 

A^ + 3 2 ^^-2 / ^\ /Ar + 6 2 2 -Ar \ 
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f N -3 1 N -2 \ 

5 2 f8 N-2 N \ (N-l y N-2 \ 
Py--ey + ^y +u^-y + -^z -jwj+v {^-^w + - - -^zj , (4.4) 

TV- 7 2 , ^ , y , N-12 



= -€z + 2wz^ —z^ + 2zu + v \^-w - I + — - — zj , (4.5) 

where u, V, w, y, and z are the MS renormahzed couphngs (see footnote 2 for a precise 
definition). 

For N < Nc ^ 21.8 + 0(e) one finds only six FPs (they belong to the models considered 
in Sec. Ill A) and none of them is stable. At N — Nc four new FPs appear: they belong to 
the 0{2)<SiO{N) models considered in Sec. Ill A [cf. cases (b) and (b')]. It is easy to check 
that none of them is stable, as it happens in three dimensions, cf. Sec. IIIB. Therefore, the 
enlargement of the symmetry to 0{2)®0{N) is never realized near four dimensions, for any 
value of TV. Stable FPs are found only for N > Ni 23.97. For Ni < N < N2 f=:i 24.15 the 
stable FP has ^; = and does not belong to any of the models considered in Sec. Ill A. For 
N > N2 the only stable FP is the decoupled chiral FP. Note that, since the zeroes of the 
/9-functions are not degenerate for all A?" > 2, no new FP can emerge within the e expansion 
at higher order. The conclusions on the stable FPs apply unchanged to the multicritical 
theory with vq = since the stable FPs always have v = 0. 

This analysis can be extended to higher orders of the e expansion, using the five-loop MS 
series. On the basis of our experience with the 0(2)(8)0(A") model we expect this analysis to 
be reliable only for large A^. Indeed, in the 0(2)(8)0(A^) theory the e expansion is not able 
to provide the correct three-dimensional scenario for the physically interesting cases N = 2 
and 3, see, e.g., Ref. [20]. Therefore, we only studied the stability of the decoupled chiral 
FP, that, according to the analysis presented in Sec. IIIC, should be stable for A^ > 6. For 
this purpose we computed the stability-boundary function A^2(e) = n2{e) + l — 24.15-|-0(e). 
Repeating the analysis of Sec. IIIC close to four dimensions, it is easy to show that n2{e) is 
determined from 

1 1 

4 - e, (4.6) 



where i'^{n2,e) and I'^pi^e) are the exponents for the chiral theory 0(2)®0(n2) and for the 
XY model, respectively. At five loops we obtain 

n2(e) = 23.1513 - 28.2072e + 18.8689e^ - 27.1434e2 + 74.2373e^ + O(e^) . (4.7) 

A Pade-Borel analysis of this series gives A^2(l) ~ 7, 4, 7, at five, four, and three loops. 
The estimates nicely oscillate and allow us to predict 4 < A^2(l) ^ 7. We also analyzed the 
inverse series l/n2(e). A Pade-Borel analysis of the series is possible only at five loops; at 
four loops all Pade approximants are defective. In this case we obtain A'2(l) ~ 8. These 
results are fully compatible with the conclusions of Sec. IIIC. 
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B. Renormalization-group flow in the 3d-MS scheme for N = 3 

In order to investigate the RG flow in the complete space of the quartic parameters 
and check for the existence of other stable FPs, we compute the five-loop series of the P 
functions associated with the five quartic couplings in the MS scheme for the general theory 
(1.3). The diagrams contributing to this calculation are a few hundreds. We handled them 
with a symbolic manipulation program, which generated the diagrams and computed the 
symmetry and group factors of each of them. We used the results of Ref. [25], where the 
primitive divergent parts of all integrals appearing in our computation are reported. The 
numerical analysis of the series is very complex and long, so that we only study the physical 
case N = 3. 

To find the stable FPs wc follow the RG flow generated by the resummed /? functions. 
For the resummation we use the Pade-Borcl method, since the large-order behavior of the 
series, needed to perform the conformal mapping, is not known. We use several different 
Pade approximants, the [4/1] and the [3/2] approximants at five loops, and the [3/1] and 
the [2/2] at four loops. We consider several initial conditions close to the Gaussian FP, 
that allow us to explore a large region in the space of the five renormalizcd couplings. 
For most of the employed approximants we found runaway trajectories, up to where the 
resummation is effective. There are however two notable exceptions. First, if we use the 
[4/1] Pade approximant, we find that the RG fiow ends at the chiral 0(2)(8)0(3) FP. This is 
in agreement with the conclusions of Sec. IIIB: this FP may be the stable FP of the model. 
Second, at four loops, by using the [3/1] Pade approximant, one observes a stable FP with 
V > 0. However, there is no indication for such a FP at five loops, even as an unstable FP. 

In conclusion, the Pade-Borel analysis of the RG fiow does not provide any evidence 
for a stable FP beside the chiral 0(2)(8)0(3) FP. Thus, the multicritical transition is either 
controlled by this FP or is of first order. 

We also considered the multicritical model with v = 0. The nonperturbative analysis of 
Sec. IIIC indicates that the decoupled collinear FP is stable for = 3. Wc also analyzed 
the full model, checking whether other FPs are present. The five-loop analysis did not find 
additional FPs, so that the multicritical behavior should be controlled by the decoupled 
collinear FP (of course, when the transition is continuous). 

V. CONCLUSIONS 

In this paper we study a general model with symmetry 0(2)C>?)[Z2©0(A^ — 1)] for > 3, 
focusing on the nature of the multicritical point. For ^ 6, the theory presents a stable 
FP that is expected to control the multicritical behavior: the decoupled FP where the fields 
(j) and become separately critical. Their critical fluctuations are controlled respectively 
by the chiral 0(2)(g)0(A^ - 1) FP and by the XY FP. For 3 < < 6 neither the decoupled 
FP nor the chiral 0(2)cg>0(A^) FP are stable. However, we have not performed a thourough 
analysis of the RG fiow, so that we cannot exclude that there exists a nontrivial stable FP 
that does not belong to any submodel we have investigated. A complete analysis has been 
performed for the physically interesting case N — 3. In this case, we find that the only 
possible stable FP is the chiral 0(2)00(3) FP. The perturbative analysis in two different 
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T T 

FIG. 3. Possible phase diagrams for easy-axis STAs in a magnetic field: (a) the multicritical 
transition is continuous; (b) the multicritical transition is of first order. Continuous thick lines 
represent first-order transitions, while the dashed lines correspond to second-order transitions. 

perturbative schemes is unable to draw a definite conclusion on the stability of this FP. 
The possible phase diagrams are reported in Fig. 3. In phase diagram (a) the transition 
at the multicritical point is continuous and controlled by the 0(2)(8>0(3) FP, thus showing 
the symmetry enlargement originally put forward in Ref. [12]; this scenario requires the 
stability of the 0(2)(8)0(3) FP and that the system is within its attraction domain. In 
the other possible cases, i.e. if the 0(2)00(3) FP is unstal)lc or if the system is outside its 
attraction domain, we should observe phase diagram (b), in which the multicritical transition 
is of first order; first-order transitions are also expected along the lines separating phases 
1-4, 1-2, and 1-3, close to the multicritical point. 

It is interesting to compare this scenario with the experimental results [1]. The behavior 
observed in experiments at the multicritical point in always compatible with a second-order 
transition. Moreover, the experimental estimates for the critical exponents (see the results 
for CsNiCla and CsMnIs in Ref. [1]) are reasonably close to the theoretical results for the 
0(2)®0(3) chiral universality class. Therefore, the phase diagram presented in Fig. 3(a) 
seems to be favored by experiments, even though the first-order scenario is not necessarily 
ruled out. If phase diagram (b) is the correct one, a possible explanation of the experiments 
is that the first-order transition is rather weak, so that all experiments are still probing a 
crossover region. This interpretation may be supported by the following reasoning. The 
experimental systems have a small easy-axis anisotropy and therefore, see Sec. II, they are 
approximately described by the effective theory with ~ yo ~ -2^0 ~ 0. Thus, the RG flow 
starts very close to the 0(2)(8)0(3) FP, so that one expects strong crossover effects controlled 
by the chiral 0(2)(8)0(3) theory. Moreover, the flow out of the chiral 0(2)00(3) FP should 
be very slow, since the associated crossover exponent </)4,4 is very small, 04,4 < 0.1. The first- 
order nature of the multicritical point also implies that along the lines separating phases 1-2 
and 2-3, one should observe first-order transitions, a tricritical point, and then XY behavior. 
The presence of the tricritical point might explain why at the transitions ior H — one 
observes values of and 7 that are significantly different from the XY estimates and that 
are closer to the mean-field predictions (see the results for CsMnIa reported in Ref. [43]). 
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Finally, the presence of the first-order transition might also explain the discrepancies between 
the experimental estimates of the critical exponents along the line 1-4 and the theoretical 
predictions for the 0(2)® 0(2) chiral universality class. Indeed, as discussed in Ref. [20], 
close to the first-order parameter region, 0(2) (8)0 (2) chiral systems show strong crossover 
effects, with effective exponents that may significantly differ from their asymptotic value. 

APPENDIX A: MEAN-FIELD PHASE DIAGRAM 
1. Model with vo = 

We begin by discussing the mean-field phase diagram for vq — 0, which corresponds to 
the model with larger symmetry [0{2)<SiO{N — 1)]©0(2). For the discussion it is useful to 
introduce new couplings 

gi^Uo + Wo, 
5-2 = Mo + 2/0, 

§3 = 2uo + 2wo - zq. (A1) 

In terms of gi, g2, and g^, the stability conditions for the quartic potential are particularly 
simple. We have gi > 0, g2 > 0, g^ > 0, and 

uo > - a/^2W2 for < g^K 2gi, 

Uq > -^Jgxgi for g-i > 2gi. 

In order to solve the mean-field equations it is useful to use the symmetry in order to obtain 
a simple parametrization of the fields. Using the 0(A'" — 1) invariance for and the 0(2) 
invariance for we can write 

0,1 = (a, 0,...), = (6,c,0...), Vi = (c?,0). (A2) 

Note that we have not used the additional 0(2) symmetry transformations applied to the 
4> fields. If a 7^ and \? \ (? ^ 0, we can perform 0(2)®0(A'" — 1) transformations to set 

The mean-field equations are easily solved and we obtain six classes of solutions (we 
report one representative for each class; other solutions in each class are obtained by applying 
the symmetry transformations): 

(a) a — h — c — d — with energy H — 0. This is a minimum only if ri > and r2 > 0. 

(b) (pl = —6ri/gi, 02 = 0, ip"^ = 0, with energy H = —\r1/gi. 

(c) — — —Gri/g^, b — d — 0, with energy H — —Srl/g^. 

(d) = -6r2/g'2, a = 6 = c = 0, with energy H = -|rf/cj'2- 

(e) = -6{g2ri - r2Uo)/D, d^ = -6{gir2 - riUo)/D, = = 0, with H = -^{g2rl + 
girl - 2rir2Uo)/D, with D = gig2 - ul. 
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(f) a?^c^ = -6{g2n-r2Uo)/D2, <f = -6{g3r2-2nuo)/ D2, with H = -|(2^2r? + y3ri - 
irir2Uo)/D2, with D2 = 9293 - 2ul, 



where ri = and r2 — r^p. In order to determine the phase diagram, for each value of the 
couphngs and of the ratio ri/r2, we must determine which solution has the lowest energy. 
The results are the following: 

(1) < 93 < 291 and < 9392/2. The critical point is tetracritical. For ri > and 
r2 > the system is paramagnetic. Then, proceeding anticlockwise in the plane 
(ri,r2) we find: a transition line for ri = and r2 > 0; phase (c); a transition line 
for r2 = ^u^rxj and ri < 0; phase (f); a transition line for r2 = 92^1/ and r2 < 0; 
phase (d); a transition line for r2 = and ri > 0. All transitions are continuous. In the 
presence of fluctuations transitions (a)/(d) and (c)/(f) belong to the XY universality 
class, while transitions (a)/(c) and (d)/(f) are chiral transitions, i.e. correspond to the 
symmetry breaking 0{2)®0{N - 1)^0{2)®0{N - 3). 

(2) < < 291 and > 5'3g'2/2. The critical point is bicritical. For ri > and r2 > the 
system is paramagnetic. Then, proceeding anticlockwise in the plane (ri,r2) we find: 
a transition line for ri = and r2 > 0; phase (c); a transition line for r2 = 

and Ti < 0; phase (d): a transition line for r2 = and ri > 0. Transitions (a)/(c) 
and (a) / (d) are continuous: transition (a) / (c) is a chiral one while transition (a) / (d) 
belongs to the XY universality class. Transition (c) / (d) is of first order. 

(3) 513 > 291 and < 9192- The critical point is tetracritical. For ri > and r2 > the 
system is paramagnetic. Then, proceeding anticlockwise in the plane (ri,r2) we find: 
a transition line for ri — and r2 > 0; phase (b); a transition line for r2 = uori/91 
and ri < 0; phase (e); a transition line for r2 = 92Ti/uq and r2 < 0; phase (d); a 
transition line for r2 = and ri > 0. All transitions are continuous. In the presence 
of fiuctuations transitions (a) / (d) and (b) / (e) belong to the XY universality class, 
while transitions (a)/(b) and (d)/(e) are coUinear transitions, i.e. correspond to the 
symmetry breaking 0(2)®0(A^ - 1)^ Z2®0{N - 2). 

(4) 5(3 > 291 and Uq > 9192- The critical point is bicritical. For ri > and r2 > the 
system is paramagnetic. Then, proceeding anticlockwise in the plane (ri,r2) we find: 
a transition line for ri — and r2 > 0; phase (b); a transition line for r2 = ri \/ 921 9\ 
and T\ < 0; phase (d); a transition line for r2 = and ri > 0. Transitions (a) /(b) and 
(a) / (d) are continuous: transition (a) / (b) is a collinear one while transition (a) / (d) 
belongs to the XY universality class. Transition (b) / (d) is of first order. 

2. Model with t-o > 

We shall now focus on the case vq 7^ 0. We shall only consider the case vo > 0, since 
the case vq < can be recovered by using the symmetry (3.2). As before, we introduce new 
couplings 
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gi^Uo + Wo, 

53 = 21^0 + 2wo -Vq- Zq, (A3) 

that are all required to be positive by the stability of the quartic potential and are invariant 
under the symmetry (3.2). We also found additional necessary stability conditions. First, 
we must have 

uo > -y/gm, (A4) 

where only the first condition is relevant for vq > 0. In order to write down the second 
condition let us define 

Ri = 7:9293 - 



2'''' 4 2^1-^3' 



2(2(71 -(73) 

R3 = (2fi'i5'2 - 9293 - 25-1^0 + 93U0 + 9m - UqVo) 

X (25(15(2 - g293 - 251^0 + 53^0 + 91^0 - 93V0 + uqVq - v^), (A5) 

and the domain Q in the coupling space 

Q = {{uo, Vq, Qi, g2, 93) ■ 93 < Sg'i, 5-2 > -^2, -^3 > O}. (A6) 

If {uq,Vo, gi, g2, 93) G Q, then the couplings must satisfy 

> Y - VRi, (A7) 

We have not been able to prove that these conditions are sufficient for the stabihty of the 
quartic potential. However, since they go over to the stability conditions for f o = (in this 
case we have proved they are sufficient), we believe they are enough for the stability of the 
quartic potential. 

We parametrize the (pai fields as before, while for tjji we must keep both components, i.e. 
we set ipi = {d, e), because of the reduced symmetry of the model. Setting as before ri = 
and r2 — r^, the mean-field equations are 

an + ^ [uo{(j)'^ + V^2) - vo{c^ + e^) + Wo0^ - z^c^] + ^vohde = 0, (A8) 

hvi + ^ \uo{(t)'^ + i)'^) - vod^ + wo0^1 + Ivoade = 0, (A9) 



cri + ^ [uo{(f)'^ + ij^) - vo{a^ + d^) + wo(/>^ - z^a^] = 0, (AlO) 
dr2 + i [uo{(f' + i^^) - Vo{b^ + c^) + y^ij^] + Ivoabe = 0, (All) 

D 

e 1 
er2 + g [uo{(p^ + - voa^ + VqiP^] + -v^ahd = 0. (A12) 

The solutions are: 
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(a) a — b — c — d — e — 0, with energy H — 0. 

(b) (f)l = -Qri/gi, 02 = 0, = q, with energy H = -^rl/gi. 

(c) = = —Qri/g3, b = d = e = 0, with energy H = —Srl/g^. 

(d) = -6r2/5'2, a = 6 = c = 0, with energy if = -|rf/g'2- 

(el) 6^ = -6(5f2?^i - T2«o)/-D, = -6(5fir2 - riUo)/D, c? ^ (? ^ ^ ^ 0, with H = 
-|(5'2r? + S'lrf - 2Tir2UQ)/D, with D = gig2 - ul- 

(e2) 6^ = -6(5(2ri - r2(Mo - ^o))/A, c?^ = -6(5'ir2 - ri{uo - vo))/D^, = = = 0, 
with H = -\{g2rl + 5'i^2 " 2rir2(Mo - wo))/^3, with = 5(15(2 - {uq - vof. 

(f) This sohition is too long to be reported. Both vectors are nonvanishing, and we can 
take a? = b'^ + (? , b^ = 0{vq). Moreover, d"^ = and de/{ab) = (53 — 2gi)/vo. 

The solutions follow the labelling used for vq = 0. Note that for vq > there are two 
different solutions corresponding the solution (e) found before. The derivation of these 
solutions is straightforward, except for case (f). To derive (f), assume that a, b, c, d, and e 
are nonvanishing. Then, Eqs. (A9) and (AlO) imply 

de = ab^-^^. (A13) 

Substitute this relation in Eqs. (A8) and (A9). Analogously, one can use this relation to 
express ab in terms of de in Eqs. (All) and (A12). The five mean-field equations become 
linear in a^, 6^, c^, d^, and e^. However, only three of them are independent. To completely 
solve the problem, we consider Eqs. (A8), (A9), and (All) that are independent, the equation 
= 62+c2 that fixes the 0(2) invariance, and the relation between rf^e^ and a^6^ that follows 
from Eq. (A13). This system of equations has a unique solution for a^, 6^, c^, d^, and e^. 

The analysis of the phase diagram of this model is extremely complex, mainly due to the 
cumbersome expressions for solution (f). We have used analytic and numerical methods to 
sort out the different possibihties. We find: 

(a) We have a pentacritical point for 251 — 5(3 > 0, g2 > R2, and Vq/2 — \/Ri < Uq < 
Vq/2 + \fRi. For ri > and r2 > the system is paramagnetic. Then, proceeding 
anticlockwise in the plane (ri, r2) we find: a transition line for ri = and r2 > 0; phase 
(c); a transition line for r\ — g^r2/{2uQ — Vq) and ri < 0; phase (f); a transition line; 
phase (e2); a transition line for ri = {uq — Vo)r2/g2 and r2 < 0; phase (d); a transition 
line for r2 = and ri > 0. All transitions are continuous except that between phases 
(f ) and (e2) that is of first order. In the presence of fiuctuations transitions (a) / (d) 
and (c)/(f) belong to the XY universality class, transition (d)/(e2) belongs to the 
0{N — 1) vector universality class, while transition (a)/(c) is a chiral transition, i.e. it 
corresponds to the symmetry breaking 0{2)^0{N — 1)— >0(2)(8)0(iV — 3). For vq ^ 
the width of phase (e2) goes to zero and we obtain case (1) considered in the previous 
Section. 
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(b) We have a tetracritical point in three cases: 2gi — > 0, g2 > R2, uq < vo/2 — a/Ri 
up to the stabihty boundary; 2gi — g^ > 0, g2 > R2, Vo/2 + \/Ri < uq < vq-\- ^J~g2g^J2\ 
2(71 — fi's > 0, < (72 < i?2, ""o < "i^o + a/ 5'2fi'3/2 up to the stabihty boundary. For ri > 
and r2 > the system is paramagnetic. Then, proceeding anticlockwise in the plane 
(ri,r2) we find: a transition line for ri = and r2 > 0; phase (c); a transition line 
with ri < 0; phase (e2); a transition line with r\ — {uq — vo)r2/g2 and r2 < 0; phase 
(d); a transition line for r2 = and ri > 0. Transitions (a)/(c), (c2)/(d), (d)/(a) are 
continuous, while transition (c) / (e2) is of first order. In the presence of fiuctuations 
transition (a)/(d) belongs to the XY universality class, transition (d)/(e2) belongs to 
the 0{N — 1) vector universality class, while transition (a)/(c) is a chiral transition, 
i.e. it corresponds to the symmetry breaking 0{2)<SiO{N — 1)— >0(2)(8)0(A'" — 3). Such 
a case does not exist for vq — 0. 

(c) We have a bicritical point for 2gi — g^ > and uq > vq + \fg2gzl2. For ri > 
and r2 > the system is paramagnetic. Then, proceeding anticlockwise in the plane 
(ri,r2) wc find: a transition line for ri = and r2 > 0; phase (c); a transition line 
for r\ = r2{gz/2g2Y^'^] phase (d); a transition line for r2 = and ri > 0. Transition 
(c) / (d) is of first order, transition (a) / (c) is a chiral transition, while transition (a) / (d) 
belongs to the XY universality class. 

(d) We have a tetracritical point for 2gi — gs < and vq — ^gig2 < Uq < Vq + y/gig2- 
For ri > and r2 > the system is paramagnetic. Then, proceeding anticlockwise 
in the plane (ri,r2) we find: a transition line for ri — and r2 > 0; phase (b); a 
transition line for ri = gir2/{uo — Vq) and ri < 0; phase (e2); ri = {uq — Vo)r2/g2 
and r2 < 0; phase (d); a transition line for r2 = and ri > 0. All transitions are 
continuous. Transition (a)/(b) is a collinear transition with Vq < 0, i.e. it corresponds 
to the symmetry breaking 0(2)(8)0(iV — 1)— > Z2® 0(A^ — 2), transitions (b)/(e2) and 
(a)/(d) are XY transitions, while transition (d)/(e2) belongs to the 0(A'" — 1) vector 
universality class. 

(e) We have a bicritical point for 2gi — gs < and uq > vq + sjg\g2- For ri > and 
r2 > the system is paramagnetic. Then, proceeding anticlockwise in the plane 
(ri,r2) we find: a transition line for ri = and T2 > 0; phase (b); a transition line for 
ri = r2{gi/g2Y^^', phase (d); a transition line for r2 = and ri > 0. Transition (a)/(b) 
is a collinear transition with vq < 0, transition (b) / (d) is of first order, and transition 
(a) / (d) belongs to the XY universality class. 

APPENDIX B: CLASSIFICATION OF THE PERTURBATIONS OF AN 

0{M)^0{N) THEORY 

We consider an 0(M)(g)0(A^)-symmetric theory with M,N > 2 and the Hamiltonian 

nMN= + r$^; 

ai 
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ai a,b i i i 

where (c^ = 1) •••^ and i = 1, ...M) are N x M matrix variables. We wish now to classify 
the quadratic and quartic operators that break the 0{N) symmetry but preserve the 0(M) 
symmetry. The classification of all the quadratic operators that break both symmetries is 
presented in App. C. 

At the quadratic level, there is only one operator, the spin-2 operator 

= $a • $6 - (B2) 

where (as in all this Section) the scalar product is in the 0(M) space and $^ = $a ■ $a- 
At the quartic level, there are three operators, one spin-4 operator, two spin-2 operators 
and one operator associated with a nontrivial 0(iV) representation. The spin-4 operator is 
given by 

(^it2 = (*a-*;>)(*c-*d) + 2perm. 

~7VT4*'^'^"^^^' • ^d) + 5 perm.) 
~7V + 4^'^"^ ■ • + 5 perm.) 



e 



+ (jV + 2KiV + 4) E(f' ■ */)'l('5<*'5o^ + 2 P™.), (B3) 

where the permutations are such to make O^^cd symmetric and traceless. The spin-2 oper- 
ators are given by 

.^(4,2,1) _ ^2^/^(2,2) 
'^ab — ^ '^ab ' 

OT"^ = • *e)($. ■ $e) - ^Sab[J2{<l>e ■ $/)^]. (B4) 

e ef 

They are always independent except for N = 2. In this case O^^l^^'^'' = 0^''^'^\ 
The remaining operator is 

= i^a ■ *c)(*6 • *d) - (*a • *d)(*6 " *c) 



__0a^. „(4,2,1) _ ^/^(4,2,2)x _ /a^(4,2,1) _ /o(4.2,2)\ 
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Such a quantity satisfies the properties Oifi = = ^25 = 

X^^i C^mcA = 0- It belongs to a nontrivial representation of the 0{N) group, the one 
associated with the Young tableau that has the shape of a 2x2 box. (Note that in our 
terminology we call spin-A; representation the representation associated with the Young 
tableau that has the shape of a row of length k). This representation occurs only for TV > 4. 
Indeed, this operator is not defined for = 2; for A?" = 3 we have '}2abcd ^^"'^^^^'''^tbld ^ 0' 
which shows that the operator does not occur. The operator C^^^^ does not mix with any 
other (quadratic or quartic) operator. 

These operators control the symmetry breaking 

0(M) ® 0{N) 0{M) ® [0{k) ® 0{N - k)] . (B6) 

The corresponding multicritical Hamiltonian is 

where = ELi Of^\ ^(^•^) = ELi C^^^, V^^^'^'^^ = ELi 0'^^^'\ V^^'^'^) = 

Ea=i V^^^-"^ = Ylh=i C'Safe- For A; = 1 or = - 1 , V^^'"-^ vanishes. Indeed, if A: = 1 

^(4,r) ^ Qi^-:^ ^i^g^^ vanishes because of the antisymmetry of the indices. For k — N — 1, 
using the fact that O^^'d traceless, we have V^'^''^^ = — Ea=i ^^aNaN — ^n'nnn — 0- 

Exphcitly, if $am ^ {4>Am, 'ipam) with A = 1,. . .k, a = 1,. . .N - k, m = 1,. . .M, we 
have 

- (^^ + 2)(iV + 4) {-^^ - '^^^ - ' + + ' • ^-)^] 

-A;(A: + 2)10;^ + 2 5^(t/>« • + 2(A: + 2)(iV - + 2)[0V + 2 $^(0a • ^a)' 

AB a/3 

-MA; - - E(^« • ^/^)'] + - 1)(^ - k - 1)[0V - J2(^A ■ V'a)1 

a/3 Aa 

It is trivial to check that y^^-'') vanishes for A; = 1 or /c = A?" — 1. 

In conclusion, at the quartic level: for A?" = 2 there are two breaking operators, V^^'^'^^ 
and for A^ = 3 and for A^ > 4, A; = 1 or A; = A^ — 1, there are three breaking operators. 
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y (4,2,1)^ yW)^ y(4,4). for AT > 4 and 2 < A; < - 2, all four operators should be 
considered. 

These operators are also relevant in other cases. If we consider the breaking 

0{M) ® 0{N) 0{M) Cn, (BIO) 
where Cat is the cubic group in an iV-dimensional space, we obtain the Hamiltonian 

N 
o=l 

For iV = 4 we may consider the breaking 

0{M)®0(A)^0(M)®SO{A). (B12) 
The corresponding Hamiltonian is 

N 

HMN + fJ2 ^"'^'^Sci- (B13) 

abcd=l 

Note, finally, that reduction to smaller symmetry groups does not require the consideration 
of additional operators, although there may be additional terms in the Hamiltonian. For 
instance, the breaking 

OiM)®0{N)^0{M)^{Ck(BCN-k) (B14) 
is obtained by considering 

k N 

n^,k + hYl ^^^al + /6 E ■ (B15) 

a=l b=k+l 

The operators proportional to /i, /5, and fe are of course degenerate at the 0(M)(g)0(A^) 
FP. 



APPENDIX C: RG DIMENSIONS OF THE QUADRATIC PERTURBATIONS AT 

THE 0(2)®0(A^) FIXED POINTS 

In this appendix we consider the three-dimensional 0{2)®0{N) invariant theory, cfr. 
Eq. (Bl) with M — 2, and compute the RG dimensions of all quadratic operators breaking 
the 0(2)®0(A^) symmetry at the 0(2)®0(A^) FPs. 

For generic M, the quadratic operators breaking the 0{M)<SiO{N) symmetry are exphc- 
itly given by 
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(CI) 

^ij Yl "^-^k^bk + -MNSabSij ^ck^ck, (C2) 
k ck 

(C3) 
(C4) 

where is a real field with a = 1, . . . ,N and i = 1, . . . , M. These operators have a simple 
group-theory intepretation, that allows us to check that the list is exhaustive. The operator 
^alj transforms as a spin-1 operator under both 0{N) and 0(Af), O^^,^- transforms as a 

spin-2 operator under both groups, while O^^-* and O^f transform as a scalar under one group 
and as a spin-2 operator under the second group. For M — 2 these operators correspond to 
those reported in Ref. [14]. 

For M = 2, we computed the RG dimensions t/j of the above-reported quadratic operators 
to six loops in the MZM scheme and to five loops in the 3(i-MS scheme. We used a symbolic 
program to generate diagrams and group factors and the compilations of Feynman integrals 
of Refs. [39,25]. We do not report the series that are available on request. The results of 
the analyses, using the resummation methods outlined in Refs. [40,19,20], are reported in 
Table III for several values of A^. In the case = 6 we do not report results for the chiral 
FP in the MZM scheme, since in this scheme there is little evidence for the existence of a 
FR The RG dimensions of the operator O^^^ related to the chiral degrees of freedom at the 
chiral FP have already been computed by exploiting the same FT methods in Refs. [44,20]; 
we report them here for the sake of completeness. We also mention that these exponents 
have been computed to order 1/N in Ref. [14], while Ref. [45] reports a calculation of 
the RG dimension of 0^^\ 

For N = 2 the RG dimensions at the collinear FP can be related to the RG dimensions 
of operators in the XY model. Indeed, the 0(2)®0(2) coUinear FP is equivalent to an XY 
FP. The mapping is the following. One defines two fields Oj and 6^, i — 1,2, and considers 
[14] 

011 = (oi - b2)/V2, 

022 = (Ol + b2)/V2, 

012 = {bi - a2)/\/2, 

021 = {bi + a2)/V2. (C5) 

At the collinear FP, fields a and b represent two independent XY fields. Using this mapping 
it is easy to show that: O^^j ~ -|- 6^; Of^j is the sum of — a^, h\ — h\, aia2, and 6162; 

(o^ O'lf) is the sum of terms of the form aibj. Thus, yi = yt^xY = I/2 = ^t.xy, 

where |/t,xy is the RG dimension of the spin-2 quadratic operator in the XY model, and 
1/3 — Hi — 2y/i,xY — 3, where yh,xY = (5 — ^xy)/2 is the RG dimension of the field in the 
XY model. Note also that the scalar operator 0^ becomes a"^ + b"^, as obviously expected. 
Estimates of i^xy, ?7xy, and yT,XY can be found in Refs. [38,22,18,42]. 



0% = -2 i^a^bi + ^a,^b^) ' j^S^b ^ 

c 

^ab = ^'^k^bk - ^Sab Y ^ck^ck, 

k ck 

c ck 
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TABLE III. RG dimensions of the quadratic operators breaking the symmetry 0{2) ® 0{N) 
at the chiral (ch) and cohinear (cl) FPs. We report the estimates obtained by analyzing the MS 
(five loops) and MZM (six loops) expansions. The errors include the spread of the considered 
approximants and the uncertainty on the location of the FP. The results for A?^ = 2 at the collinear 
FP have been obtained by using the mapping with the XY model and the numerical results reported 
in Refs. [38,18]. 



FP.AT 




yi 




2/2 




ys 




2/4 




MS 


MZM 


MS 


MZM 


MS 


MZM 


MS 


MZM 


ch,2 


2.37(15) 


2.54(12) 


2.00(15) 


2.07(7) 


1.34(15) 


1.25(4) 


1.34(15) 


1.25(4) 


ch,3 


2.25(12) 


2.35(13) 


1.96(11) 


1.99(4) 


1.54(8) 


1.49(3) 


1.21(9) 


1.09(5) 


ch,4 


2.17(10) 


2.29(8) 


1.94(10) 


2.04(20) 


1.65(3) 


1.64(5) 


1.17(8) 


1.06(7) 


ch,5 


2.05(10) 


2.20(7) 


1.93(10) 


1.98(20) 


1.72(4) 


1.72(8) 


1.15(7) 


1.02(8) 


ch,6 


2.03(7) 




1.90(10) 




1.76(4) 




1.13(9) 




ch,8 


2.02(2) 


2.03(4) 


1.92(2) 


1.93(4) 


1.81(2) 


1.79(1) 


1.13(8) 


1.13(4) 


ch,16 


2.001(5) 


2.00(1) 


1.948(5) 


1.95(1) 


1.897(7) 


1.885(5) 


1.08(2) 


1.07(1) 


chjOO 




2 




2 




2 




1 


cl,2 




= 1.489(6) 




= 1.766(6) 


21/;,, XY 


- 3 = 1.9620(8) 


2!//t,XY 


- 3 = 1.9620(8) 


cl,3 


1.2(1) 


1.15(10) 


1.75(5) 


1.75(10) 


2.0(1) 


2.0(2) 


2.1(2) 


2.05(15) 


cl,4 


1.1(1) 


1.10(15) 


1.65(10) 


1.66(5) 


1.90(15) 


1.75(10) 


2.0(3) 


2.05(15) 
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